Evolution of Hydraulic Conductivity of Unsaturated Compacted Na-Bentonite under Confined Condition—Including the Microstructure Effects

Compacted bentonite is envisaged as engineering buffer/backfill material in geological disposal for high-level radioactive waste. In particular, Na-bentonite is characterised by lower hydraulic conductivity and higher swelling competence and cation exchange capacity, compared with other clays. A solid understanding of the hydraulic behaviour of compacted bentonite remains challenging because of the microstructure expansion of the pore system over the confined wetting path. This work proposed a novel theoretical method of pore system evolution of compacted bentonite based on its stacked microstructure, including the dynamic transfer from micro to macro porosity. Furthermore, the Kozeny–Carman equation was revised to evaluate the saturated hydraulic conductivity of compacted bentonite, taking into account microstructure effects on key hydraulic parameters such as porosity, specific surface area and tortuosity. The results show that the prediction of the revised Kozeny–Carman model falls within the acceptable range of experimental saturated hydraulic conductivity. A new constitutive relationship of relative hydraulic conductivity was also developed by considering both the pore network evolution and suction. The proposed constitutive relationship well reveals that unsaturated hydraulic conductivity undergoes a decrease controlled by microstructure evolution before an increase dominated by dropping gradient of suction during the wetting path, leading to a U-shaped relationship. The predictive outcomes of the new constitutive relationship show an excellent match with laboratory observation of unsaturated hydraulic conductivity for GMZ and MX80 bentonite over the entire wetting path, while the traditional approach overestimates the hydraulic conductivity without consideration of the microstructure effect.


Introduction
Bentonite is a widely considered swelling clay as an engineered barrier material in geological waste disposal because it has a large percentage of smectite (50~90%), a clayey soil swelling under water intrusion. Bentonite has lower hydraulic conductivity and diffusivity, and higher expansive capacity, specific surface area, cation-exchange potential, and thermal conductivity compared to other soils [1][2][3][4][5][6]; therefore, it plays a key role in various geo-environmental engineering applications, such as containment systems, carbon dioxide storage [7,8], bioremediation [9], and stability of petroleum reservoirs [10]. Bentonite has been attracting attention in the construction of disposal repository of highlevel nuclear waste for 4 decades [11]. Sodium bentonite such as GMZ and MX80 bentonite is commercialised bentonite used more extensively than calcium bentonite in engineering because of its more outstanding expansive ability and lower hydraulic conductivity to water transport [4]. The bentonite could have side effects on personal health because it can prevent digestion and influence the absorption of electrolytes [12]. High levels of germs or heavy metals from bentonite are harmful to health [12]. Meanwhile, the wetting of bentonite clay can damage roadways and buildings without proper pretreatment [13].
In engineering applications, bentonite is normally densely compacted with very low water content, and thus it is unsaturated with a rather high suction at the beginning. Bentonite is placed between the radioactive waste and the host rock, serving as buffer/backfill material. Thereafter, bentonite consistently experiences water intrusion from the host rock. The stiffness and strength of host rock are usually too high to deform, so the wetting path of bentonite happens in a quasiconstant volume condition [14]. Having a solid understanding of hydraulic behaviour of engineering buffer (bentonite) has practical significance in a real disposal repository of HLRW that includes a confined condition and wetting path [14][15][16][17][18][19][20]. Laboratory observations found the water flow and ion diffusion are considerably influenced by the expansion and shrinkage of the microstructure of clays. The application of the conventional constitutive relationship of unsaturated hydraulic conductivity might overestimate the rate of water flow and ion transport [21,22]. Limited experimental and theoretical approaches can provide a solid demonstration of microstructure effects on unsaturated hydraulic conductivity of compacted bentonite. The traditional method for calculating the unsaturated hydraulic conductivity of soils is the product of relative hydraulic conductivity and saturated hydraulic conductivity [23][24][25][26][27][28][29][30][31]; however, both variables are strongly impacted by the swelling microstructure of compacted bentonite [23,32].
Experimental studies show that saturated hydraulic conductivity of porous media is decided by the percentage of smectite, fluid, temperature, porosity, pore pressure and pore geometry that consists of pore size distribution, tortuosity, pore throats, coordination number of pore, etc. [33][34][35][36][37][38]. Despite many efforts to describe saturated hydraulic conductivity [30,[39][40][41][42][43], a solid theoretical model of hydraulic properties of swelling clays remains a challenge. A well-developed constitutive model to estimate saturated hydraulic conductivity of coarse-grain and nonswelling soils is the Kozeny-Carman (KC) relationship [44,45]. The KC equation interprets the saturated hydraulic conductivity as a relationship of porosity (Φ), specific surface area (S A ), tortuosity (τ) and a shape factor (C s ), shown as Equation (1) below, 2 (1) where k sat is the hydraulic conductivity of saturated soils (m/s), Φ is the porosity of soils (dimensionless), C s is a dimensionless shape constant (dimensionless), S A is the specific surface area of soil (m 2 /g), γ w is unit weight of fluid (N/m 3 ), ρ d is the dry density of soil (kg/m 3 ), and η is fluid viscosity (N·s/m 2 ). The KC equation was proposed on the hypothesis of fluid through the uniform channels of a cross-section, and the laboratory observation of coarse-grained soils (e.g., sand) was in line with the results calculated by the KC equation [42,45,46]. Since the KC equation assumes no change of soil fabric and no water transport in the solid phase of particles during the wetting path, it is not valid for clayey soils [28,39,47]. The KC equation suggests that saturated hydraulic conductivity (k sat ) versus Φ 3 /(1 − Φ) 2 should be a straight line for soils. Nonetheless, experimental results failed to identify such a linear relationship in clayey soils without consideration of the microstructure effects [28].
To obtain the saturated hydraulic conductivity of swelling soils, the microstructure evolution needs to be introduced into the KC equation, mainly including the effective (external) specific surface area, effective (macro) porosity and tortuosity for water flow. Macro pores consist of interparticle and interaggregate pores for compacted clays. Diverse approaches are applied to measure the specific surface such as Grain Size Distribution Curve [47,48], despite the fact that they are not widely employed in soil mechanics and engineering hydraulics. Such approaches work for granular soils without nonswelling fine particles. Although Chapuis and Aubertin [47] presented a new method to estimate specific surface area that fit experimental data very well, it is inappropriate for clays containing smectite. The discrepancy can be ascribed to the swelling capability of clayey soils, resulting in a considerable decline of the external specific surface. In this paper, the stacked lamellar structure of bentonite is further developed based on Holmboe et al. [49] and Tournassat et al. [50], to estimate effective (external) specific surface area for water flow.
The effective (macro) porosity of bentonite is not constant because of the expansion of microstructure after water absorption. Very limited studies have been found in the theoretical description of pore system evolution of wetting bentonite. Likos and Wayllace [51] as pioneers demonstrated a simplified geometrical method to derive the macro porosity, while Sedighi and Thomas [52] proposed a novel approach to describe macro porosity by a geochemical model. To investigate the evolution of macro porosity, the interaction between water and bentonite needs to be clarified. There are three types of water states in compacted bentonite, which incorporates interlayer water located between the stacked Tetrahedral-Octahedral-Tetrahedral layers (TOT layers) of bentonite, water in diffuse double layers and free water [53]. The interlayer water is immobile, while water in the interparticle and interaggregate pores can transport freely. Water in diffuse double layers is partly restrained with much larger mobility than the interlayer water. During the wetting path, bentonite will absorb water, most of which enters interlayer pores of particles, leading to the decease of macro porosity. To predict the saturated hydraulic conductivity of compacted clays, the macro porosity needs to be applied instead of total porosity according to the realistic description of water states in the pore system of compacted clays.
Two well-known approaches employed to derive the relative hydraulic conductivity were proposed by Mualem [54] and Brooks and Corey [55] in this field. Mualem [54] described relative hydraulic conductivity as a relationship of volumetric water content as shown in Equation (2), while Brooks and Corey [55] suggested relative hydraulic conductivity as a function of suction (Equations (3)- (5)). Both of the traditional models ignore the microstructural expansion of compacted bentonite, leading to an overestimation of hydraulic conductivity and failure in explanation of the decrease of hydraulic conductivity in the early stage of wetting path (Ye et al 2009;Cui et al, 2008;Chen et al, 2020). Subsequently, van Genuchten [24] proposed a renowned water retention function that makes suction and volumetric water content convertible. In experimental observations, the hydraulic conductivity of compacted bentonite against suction shows a U-shaped curve that cannot be captured by the two traditional methods with a monotonic function; therefore, the microstructure effect needs to be taken into account to explain the observation, where k r is relative hydraulic conductivity, k sat is the saturated hydraulic conductivity, k unsat is the unsaturated hydraulic conductivity, S e is the effective saturation, θ and θ r are the actual and the residual volumetric water content, respectively, and α is a constant parameter, which is assumed to be 3.5 as an average widely adopted for soils in the literature [55,56], where µ is a dimensionless number, Ψ is suction and Ψ d is the air entry value, and λ is a fitting factor related to pore-size distribution. In this paper, compacted bentonite is applied to explore the microstructure evolution of expansive clays. This work aims to address the challenge in the overestimation of hydraulic properties of bentonite because of the pore network evolution and explore the theoretical method to describe the U-shaped hydraulic conductivity by considering both the microstructure evolution and dropping gradient of suction. The new theoretical model is proposed to describe the key hydraulic parameters with the consideration of microstructure effects, including the evolution of the specific surface area, porosity and tortuosity over wetting path. The investigation of pore evolution starts with the clay-water-chemical interaction in the stacked structure of compacted bentonite. The dynamic process of water absorption is discussed to estimate the amount of water in interlayer pores and macro pores and pore network evolution over the hydration process. The revised KC equation is proposed for predicting the saturated hydraulic conductivity of compacted bentonite by introducing the modified hydraulic parameters. The relative hydraulic conductivity is further developed to form a constitutive relationship that consists of the influence of suction and microstructure effects, followed by the comparison with experimental results.  Figure 1 shows a typical lamellar structure of one bentonite particle that consists of 10 to 350 TOT layers (Tetrahedral-Octahedral-Tetrahedral layers). The TOT layers of bentonite are stacked because of repulsive and attractive forces, forming crystal particles. The external surface area (A ext ) contributes to overall flux instead of total surface area (A tot ). A tot that is independent of water content can be measured by experimental techniques such as Ethylene Glycol Monoethyl Ether. The relationship between A ext and A tot was given as Equation (6) according to the lamellar structure of smectite [50,[57][58][59][60],

Theory and
where n is the stacked TOT layers of 1 particle with a range from 10 to 350 depending on water content and the type of bentonite sample.
structure effects, including the evolution of the specific surface area, poro osity over wetting path. The investigation of pore evolution starts with t chemical interaction in the stacked structure of compacted bentonite. The cess of water absorption is discussed to estimate the amount of water in in and macro pores and pore network evolution over the hydration process. T equation is proposed for predicting the saturated hydraulic conductivity bentonite by introducing the modified hydraulic parameters. The relative ductivity is further developed to form a constitutive relationship that consi ence of suction and microstructure effects, followed by the comparison wit results.

Pore Evolution of Compacted Bentonite
2.1.1. Specific Surface Area Figure 1 shows a typical lamellar structure of one bentonite particle 10 to 350 TOT layers (Tetrahedral-Octahedral-Tetrahedral layers). The TO tonite are stacked because of repulsive and attractive forces, forming crysta external surface area (Aext) contributes to overall flux instead of total surface that is independent of water content can be measured by experimental tech Ethylene Glycol Monoethyl Ether. The relationship between Aext and Ato Equation (6)  where n is the stacked TOT layers of 1 particle with a range from 10 to 350 water content and the type of bentonite sample.  To estimate the evolution of pore structure, the surface area (A s ) of a single TOT layer in particles (A s ) needs to be assessed. The external specific surface area is defined as the external surface of a particle divided by the mass of the particle, where m single is the average mass of per particle. The total number of TOT layers (N) in a bentonite sample is described as Equation (8).
where m smectite is the mass of smectite (kg). Therefore, the m single is expressed as below, Substituting Equation (9) to Equation (7) yields, Therefore, the combination of Equations (6) and (10) yields A s , 2.1.2. Evolution of Pore System Figure 2 presents a conceptual diagram of the pore system of compacted bentonite. The microstructure of compacted bentonite is composed of a cluster of TOT layers that constitute bentonite particles and aggregates [52,61]. In many natural porous materials, such as tight rocks or compacted swelling clays, the intraparticle and interparticle pore sizes are in the range of micro (less than 5 nm), meso (between 5 and 50 nm), and macro pores (more than 50 nm) [62]. Several different pore definitions for compacted bentonite were displayed in different publications [19,63,64] where the interlayer pores are identical to intraparticle pores and interparticle pores referred as to intra-aggregate pores. Water flow is constrained in interlayer pores because of interaction between water and the TOT surface of compacted bentonite; therefore it reaches a consensus that the water in interlayer pores is immobile [53,65]. By contrast, water can freely flow through interparticle and interaggregate pores, which contribute to long-range water transport. The interparticle and interaggregate pores are collectively referred to as macro pores, as shown in Figure 2. In laboratory observations, the number of stacked water layers between two adjacent TOT layers was summarised as one, two, three and even four with higher water content, as shown in Tables 1 and 2, in which the basal spacing is the summation of the thickness of one TOT layer and thickness of water layers between TOT layers (see the upper-right panel of Figure 2). The four water layers have not been extensively measured and accepted until now, which was found by Saiyouri et al. [66] using MX-80 bentonite. Table 1 shows In laboratory observations, the number of stacked water layers between two adjacent TOT layers was summarised as one, two, three and even four with higher water content, as shown in Tables 1 and 2, in which the basal spacing is the summation of the thickness of one TOT layer and thickness of water layers between TOT layers (see the upper-right panel of Figure 2). The four water layers have not been extensively measured and accepted until now, which was found by Saiyouri et al. [66] using MX-80 bentonite. Table 1 shows that the thickness of a single TOT layer is mostly reported as 9.5 Å and the single water layer thickness is 3 Å. The diameter of a water molecule is reported as 3 Å at standard temperature and pressure [67,68], which is identical to the thickness of a single water layer presented in Table 1. Accordingly, the maximal basal spacing is the summation of three water layers and one TOT layer, equal to 18.5 Å, which is consistent with the summary of Table 1. The asterisk (*) represents the data that were estimated by researchers rather than by their experimental observations.
The adsorption reaction rate of interlayer water is comparable with hydrodynamic adsorption, which is presented by a pseudo-first order kinetic model [73]. Therefore, the evolution of basal spacing against water content is developed based on the hydrodynamic adsorption equation, where q e and q t (kg/kg) are the mass of absorbed solute at equilibrium and at time t (h), and k a (h −1 ) is the rate constant of the pseudo-first order kinetic adsorption. As stated by the boundary conditions (t = 0, q t = 0; t = t, and q t = q t ), the integral of the hydrodynamic adsorption equation yields, If the water layers are stacked in order then, If the water content (ω, kg/kg) is assumed to have a linear relationship with time, then yields, where b is a constant correlating with water absorption rate, which is −5.0 obtained by the experimental data (see Figure 3), d TOT and d t are the thickness of a single TOT layer and the basal spacing, respectively, e is natural logarithm.
If the water content (ω, kg/kg) is assumed to have a linear relationship wi then yields, where b is a constant correlating with water absorption rate, which is −5.0 obtaine experimental data (see Figure 3), dTOT and dt are the thickness of a single TOT la the basal spacing, respectively, e is natural logarithm.  The water content can be converted into volumetric water content, where ms is the mass of solid (kg), mw is the mass of water (kg), ω is the water (kg/kg), ρw is the density of water (kg/m 3 ), Vt is the total volume of the sample, V volume of total water in the sample (m 3 ). In a bentonite particle, the thickness of interlayer water between two TOT lay can be calculated as follows, The water content can be converted into volumetric water content, where m s is the mass of solid (kg), m w is the mass of water (kg), ω is the water content (kg/kg), ρ w is the density of water (kg/m 3 ), V t is the total volume of the sample, V w is the volume of total water in the sample (m 3 ). In a bentonite particle, the thickness of interlayer water between two TOT layers (d iw ) can be calculated as follows, The volume of immobile water (namely, the volume of interlayer water, V imw ) can be estimated, Consequently, the volume of macro pores (V macro ) is derived, where V s is the volume of solid. Therefore, the macro porosity (Φ macro ) can be calculated, Accordingly, where Φ macro is the macro porosity (i.e., effective porosity), Φ micro is the micro porosity (i.e., porosity of interlayer water), ρ s d is the dry density of the smectite (kg/m 3 ), G s is the specific gravity of soil, ρ d is the dry density of the soil (kg/m 3 ).
In summary, the macro porosity only correlates with one unknown variable, i.e., volumetric water content, where ρ i is a constant that was recommended 2.8 kg/m 3 by literature [52,74].

Model for Unsaturated Hydraulic Conductivity
In the hypothesis of this work, the saturated hydraulic conductivity of compacted bentonite (k sat ) can be derived by taking into account microstructure effects, where γ is the unit weight of the fluid (N/m 3 ) and η is the fluid viscosity (N s/m 2 ) [39,86].
Shape factor (C s ) is reported as 0.2 for soils [33,42,45,87,88]. The unsaturated hydraulic conductivity (k unsat ) is the product of relative hydraulic conductivity (k r ) and saturated hydraulic conductivity (k sat ) [55], The combination of Equations (2)-(5) yields the relative hydraulic conductivity as follows, The pore geometry of nonswelling soils is considered a constant whose unsaturated hydraulic conductivity is only impacted by the change of suction. Nevertheless, in reality, the unsaturated hydraulic conductivity of bentonite will be considerably affected by changing hydraulic parameters, such as macro porosity and external specific surface area and tortuosity, caused by the microstructure evolution of the pore network. In early studies, there was a lack of a cutting-edge laboratory device to measure the unsaturated hydraulic conductivity of low-permeability clays such as bentonite. The recent observations showed unsaturated hydraulic conductivity against suction should be a U-shaped curve instead of a monotonically decreasing function [15,89,90]. Therefore, this work considers swelling effects, differing from the traditional constitutive relationship such as Equations (30) and (31) where only suction influences the hydraulic conductivity. If there are swelling effects, a reference point where hydraulic conductivity is equal to fully saturated hydraulic conductivity can be found [23], The combination of Equations (31) and (32) yields, The  [15], where the suction is 35 MPa and corresponding k r = 1 as the reference point used by Liu et al. [23].
The expansive capability results in the decrease of macro porosity. Leverett [91] proposed that the change of air entry value (Ψ d ) can be represented by a function of porosity under confined condition. Therefore, Liu et al. [23] deduced Equations (34) and (35) to establish the relationship of Ψ d and air entry value at the reference point (Ψ d,ref ) under the confined condition as below, where parameter σ is larger than one for swelling materials because the porosity may not reflect the effect of well-connected pores on hydraulic conductivity. Substituting Equations (34) and (35) into Equation (33) yields, Comparing Equation (36) with Equations (34) and (35), the revised model introduces one more factor related to porosity, i.e., (φ/φ re f ) (σ−1)(1+1/µ)−3λ which represents the expansive capability of clayey soil under constant-volume condition, and (ψ/ψ re f ,c ) −(1/µ+1+3λ) reflects the change of hydraulic conductivity resulting from declining suction.
Pham et al. [92] reported that the swelling capability of clays was an exponential function of suction under confined condition, where β is a fitting factor. The combination of Equations (36) and (37) yields the constitutive relationship of unsaturated swelling clays as below, where the experimental data reported in Cui et al. [15] support that µ equals 0.28 as a constant with various suction; β can be derived by Equation (37), while The σ and λ can be found using experimental data.

Experimental Data and Parameters of Model
In the work, the experimental materials are employed to verify the theoretical model including Na-bentonite GMZ and MX80 whose properties are listed in Table 3. External specific surface area is the first parameter estimated to derive the saturated hydraulic conductivity based on the revised Kozeny-Carman equation. The experimental data of stacked TOT layers per particle were measured by Saiyouri et al. [66]. A fit to the data was proposed by Chen et al (2020), showing a good match, n = n max (1 + (n max /n min − 1)e cψ ) (39) where n max and n min are the maximum and minimum stacked TOT layers per particle with 10 and 350, respectively, given by Saiyouri et al. [66], and c is a dimensionless fitting parameter that is equal to −0.17 for Na-bentonite according to Figure 4 and is associated with the reaction rate between water and bentonite. where nmax and nmin are the maximum and minimum stacked TOT layers per parti 10 and 350, respectively, given by Saiyouri et al. [66], and c is a dimensionless fit rameter that is equal to −0.17 for Na-bentonite according to Figure 4 and is associat the reaction rate between water and bentonite.  To derive macro porosity (effective porosity), the data from Table 3 are sub into Equation (26), and the results are shown in Figure 5 (the red line). To eval accuracy, this work compares with the theoretical method of micro porosity prop Sedighi and Thomas [52] as Equation (41), where Xhs denotes the mole fraction of hydrated smectite, which can be gaine Sedighi and Thomas [52]; Msm is the molar mass of dry smectite (kg/mol), given as kg/mol by Gailhanou et al. [94]; ζc is the number of moles of water in the interla sorption or desorption reaction, reported as 4.5 mol of water, if a maximum of tw olayers of adsorbed water occurs in the interlayer pores [52]; and υil is the molar of the interlayer water, recommended as 17.22 m 3 /mol [95]. The outcomes from the two presented methods are compatible for compacted bentonite. However, the proposed method has an advantage in feasibility compar the geochemical model proposed by Sedighi and Thomas [52] because of fewer a tions during the process of derivation. Moreover, in comparison with two wate between TOT layers in Sedighi and Thomas [52], this present method takes into three water layers that are closer to the recent experimental observation for Na-be Consequently, the external specific area can be illustrated, To derive macro porosity (effective porosity), the data from Table 3 are substituted into Equation (26), and the results are shown in Figure 5 (the red line). To evaluate its accuracy, this work compares with the theoretical method of micro porosity proposed by Sedighi and Thomas [52] as Equation (41), where X hs denotes the mole fraction of hydrated smectite, which can be gained from Sedighi and Thomas [52]; M sm is the molar mass of dry smectite (kg/mol), given as 378.787 kg/mol by Gailhanou et al. [94]; ζ c is the number of moles of water in the interlayer adsorption or desorption reaction, reported as 4.5 mol of water, if a maximum of two monolayers of adsorbed water occurs in the interlayer pores [52]; and υ il is the molar volume of the interlayer water, recommended as 17.22 m 3 /mol [95].  Figure 5 shows the ɸmicro increases with a decrease in suction and finally approa stable at 0.260, while the ɸmacro drops before stable at 0.101. In the fully saturated condi ɸmicro is 2.6 times higher than ɸmacro that only makes up 27.9% of the total porosity (0.3 The outcomes from the two presented methods are compatible for compacted MX-80 bentonite. However, the proposed method has an advantage in feasibility compared with the geochemical model proposed by Sedighi and Thomas [52] because of fewer assumptions during the process of derivation. Moreover, in comparison with two water layers between TOT layers in Sedighi and Thomas [52], this present method takes into account three water layers that are closer to the recent experimental observation for Na-bentonite. Figure 5 shows the Φ micro increases with a decrease in suction and finally approaches stable at 0.260, while the Φ macro drops before stable at 0.101. In the fully saturated condition, Φ micro is 2.6 times higher than Φ macro that only makes up 27.9% of the total porosity (0.361). As a summary, as fewer than 30% of pores are macro pores that contribute to long-term water flow, the traditional KC equation using total porosity will largely overestimate the saturated hydraulic conductivity of bentonite. Secondly, the difference between macro porosity and total porosity explains clayey soils such as bentonite with a large porosity but much lower permeability compacted with nonexpansive soils. Thirdly, the transfer from Φ macro to Φ micro over a large range of suction gives rise to a decrease of unsaturated hydraulic conductivity, and as stated above, Φ micro makes no contribution to long-range water transport.

Saturated Hydraulic Conductivity
Since only macro pores contribute to water flow in compacted bentonite, the total porosity in the Kozeny-Carman equation is replaced by the macro parameter that includes external specific surface area, macro porosity and tortuosity. These macro parameters change with the swelling microstructure of compacted bentonite in the wetting path: The variables or constants from Table 3 are substituted to Equation (42), which yields that k sat of GMZ equals 1.01 × 10 −13 (m/s). Based on the macro porosity derived from Sedighi and Thomas [52], the k sat of GMZ is 1.88 × 10 −13 (m/s), whereas the experimental result of k sat is 1.18 × 10 −13 [14]. The calculated results fall within the range from 1/3 to 3 times of experimental values, which are regarded as acceptable outcomes because of the uncertainty of experimental hydraulic conductivity [28,33,86,87,96]. In summary, both methods achieve reliable results of estimated saturated hydraulic conductivity that fall within the acceptable experimental range, whereas the predictive results of the present method are closer to experimentally measured values.

Unsaturated Hydraulic Conductivity
In this section, the relative hydraulic conductivity (k r * ) of Na-bentonite including GMZ bentonite and MX80 bentonite is discussed in order to predict the unsaturated hydraulic conductivity. The reference point where k unsat = k sat has not been found for GMZ bentonite by experiments (Ye et al., 2009); however, it can be estimated that when suction is approaching 100 MPa (Ψ ref ), k unsat equals k sat , based on the trend of experimental hydraulic conductivity. Since the suction can be presented by volumetric water content according to van Genuchten [24], the macro porosity is derived using the volumetric water content through Equation (26). As van Genuchten [24] reported, the θ r of clayey soil can be regarded as zero, and the θ s of GMZ is 0.425 (m 3 /m 3 ) derived from experimental data of Ye et al. [14]. The fitting parameters "a" and "m" can be acquired from Figure 6: where β is estimated as 0.537 from Figure 7. Since the Ψref (100 MPa) is obtained by extending the trend of the experi it could be slightly different from the real value of Ψref; therefore, a small co (C × Ψ/Ψref) is introduced into Equation (38), yielding Equation (46), whic where β is estimated as 0.537 from Figure 7.
where β is estimated as 0.537 from Figure 7. Since the Ψref (100 MPa) is obtained by extending the trend of the exper it could be slightly different from the real value of Ψref; therefore, a small co (C × Ψ/Ψref) is introduced into Equation (38), yielding Equation (46), whic  (46), which reduces the deviation caused by the assumed Ψ ref . If Ψ ref can be measured accurately, the small correction term can be considered as zero, which has been proved by Liu et al. [23]. In summary, the C can be understood as a factor to express the relative error between the estimated Ψ ref and the actual Ψ ref: When the experimental data from Ye et al. [14] are substituted into the constitutive model, the λ, σ and C are equal to −1.29, 1.318, −0.64 as shown in Equation (48). From the predicted curve of Figure 8, the degree of correlation between experimental data and prediction is 99.0%, which shows the high accuracy of the theoretical model on k unsat . This model is also verified to predict k unsat using the MX80 bentonite measured by Wang et al. [90]. As shown in Figure 9, the results confirm that the present model can precisely predict the k unsat of MX80 with a degree of correlation as 99.3%.
The proposed constitutive relationship incorporates the effect of microstructural expansion on the evolution of unsaturated hydraulic conductivity. The outcome of the proposed model is in line with the experimental measurement for both bentonites. The traditional approach (Equations (30) and (31)), which ignores the microstructure effects, yields much higher prediction and is not able to explain the decrease of hydraulic conductivity with declining suction. This developed constitutive relationship can well describe water flow in swelling bentonite that has the U-shaped hydraulic conductivity. This proposed relationship reveals that macro pores (Φ macro ) continue to convert into micro pores (Φ micro ), and the average pore size of bentonite decreases over wetting path (dropping suction) of compacted bentonite, resulting in a drop of hydraulic conductivity. Subsequently, the hydraulic conductivity is dominated by the gradient of suction, showing an increasing trend in the low suction range. As a summary, the predictive outcomes of the new constitutive relationship show an excellent match with laboratory observation of unsaturated hydraulic conductivity, while the traditional approach overestimates the hydraulic conductivity without consideration of microstructure effect.

Summary and Conclusions
This work aims to develop a new constitutive relationship for unsaturated hy conductivity of compacted bentonite, including the microstructure effect. The micr ture effect is represented by the change of accessible porosity (i.e., macro poros water, further resulting in a varying tortuosity, specific surface area and relative hy conductivity. The modified Kozeny-Carman equation was developed for satura draulic conductivity incorporating the microstructure evolution, yielding a good ment with experimental results. Furthermore, a new constitutive relationship for hydraulic conductivity was developed, including both effects of capillary pressu microstructural expansions during wetting path. The prediction of the new cons relationship is in excellent accord with the laboratory observation for GMZ and bentonite over the entire wetting path, explaining the initial decrease of hydrau ductivity with declining suction (i.e., U-shaped curve) for compacted bentonite. T ditional model ignores microstructure expansion, leading to an overestimation of h lic conductivity and failure in explanation of decrease of hydraulic conductivity

Summary and Conclusions
This work aims to develop a new constitutive relationship for unsaturated hydraulic conductivity of compacted bentonite, including the microstructure effect. The microstructure effect is represented by the change of accessible porosity (i.e., macro porosity) for water, further resulting in a varying tortuosity, specific surface area and relative hydraulic conductivity. The modified Kozeny-Carman equation was developed for saturated hydraulic conductivity incorporating the microstructure evolution, yielding a good agreement with experimental results. Furthermore, a new constitutive relationship for relative hydraulic conductivity was developed, including both effects of capillary pressure and microstructural expansions during wetting path. The prediction of the new constitutive relationship is in excellent accord with the laboratory observation for GMZ and MX80 bentonite over the entire wetting path, explaining the initial decrease of hydraulic conductivity with declining suction (i.e., U-shaped curve) for compacted bentonite. The traditional model ignores microstructure expansion, leading to an overestimation of hydraulic conductivity and failure in explanation of decrease of hydraulic conductivity in the early stage of wetting path. This work helps improve the understanding of accurate prediction of flow process in the environmental application of swelling clays, including containment systems, especially in the geological disposal of high-level nuclear waste.

Informed Consent Statement: Not applicable.
Data Availability Statement: All data, models, or code that support the findings of this study are available from the corresponding author.
Acknowledgments: Tian Chen appreciates the support of the President's Doctoral Scholarship and the supervisory team from the University of Manchester.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations k
Hydraulic permeability (m/s) k r Relative hydraulic conductivity k r * Proposed relative hydraulic conductivity of compacted bentonite k sat The saturated hydraulic conductivity (m/s) k unsat The unsaturated hydraulic conductivity (m/s) S e The effective saturation a, m Curve-fitting parameters by van Genuchten [24] n max , n min The maximum and minimum number of TOT layers of bentonite particle respectively